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We examine neutrino oscillations in a two Higgs doublet model (2HDM) in which the second 
doublet couples only to the third generation right-handed up-fermions, i.e., to ta and to -/V3 which is 
the heaviest right-handed Majorana neutrino. The inherently large tan (3 of this model can naturally 
account for the large top mass and, based on a quark-lepton similarity ansatz, when embedded into 
a seesaw mechanism it can also account for the observed neutrino masses and mixing angles giving a 
very small Ova- —0.017 ^ #13 0.021 at 99% CL, and a very restrictive prediction for the atmospheric 
mixing angle: 42.9° ^ a tm ~ 45.2° at 99% CL. The large value of tan [3 also sets the mass scale of 
the heaviest right-handed Majorana neutrino N3 and triggers successful leptogenesis through a CP- 
asymmetry in the decays of the Ni (lightest right-handed Majorana) which is tan 2 (3 enhanced 
compared to the CP-asymmetry obtained in models for leptogenesis with one Higgs doublet or in 
the MSSM. This enhancement allows us to relax the lower bound on Mn ± and consequently also 
the lower bound on the reheating temperature of the early universe. 
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A monumental discovery of the past decade is that 
neutrinos have mass! This discovery is bound to have a 
significant impact on our understanding of the universe. 
It also presents us with a pressing challenge: unraveling 
of the mixing matrix in the lepton sector. When we reca- 
pitulate the analogous quark case we can easily grasp the 
importance and the difficulties of this new task. Indeed, 
an obvious but nonetheless crucial question is the rela- 
tion of the mixing matrix in the lepton sector to that in 
the quark sector. Intensive experimental effort is under- 
way and much more is being planned, round the globe, 
to address these issues of vital importance. 

In this work we suggest that these impressive findings 
in the lepton sector are closely related to another re- 
markable finding of the 90's in the quark sector, namely 
that the top quark is enormously heavy compared to 
all the other quarks. Specifically, we will propose a 
simple, grounds-up approach which treats the 3rd gen- 
eration neutrino in a completely analogous manner to 
the top quark. This reasoning leads us to a concrete 
framework for neutrino masses and mixings with con- 
siderable predictive power. For example, we find that 
the crucial mixing angle #13, which is an important tar- 
get of many neutrino experiments, is quite constrained, 
-0.017 < 0i3 £ 0.021 at 99% CL, in our picture. 

Needless to say, these observations of neutrino oscilla- 
tions pj, confirming the existence of massive neutrinos, 
implies that the Standard Model (SM) can no longer be 
regarded as a minimal theory that explains all observed 
phenomena in Particle Physics. Bearing this in mind, 
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our objective here is to present a simple extension of 
the SM, which will serve as a low energy effective theory 
which captures the dominant features of phenomenology 
and models some underlying complicated dynamics of 
electroweak (EW) symmetry breaking at distances much 
shorter than the EW scale. 

Recall the allowed 3er ranges on the atmospheric and 
solar neutrino square mass differences and mixing angles 
(including the CHOOZ v e disappearance experiment) 0] : 

7.3 < Am 2 sol ■ 10 5 eV 2 < 9.3 , 0.28 < tan 2 sol < 0.6 , 
1.6 < Am 2 tm ■ 10 3 eV 2 < 3.6 , 0.5 < tan 2 atm < 2.1 , 
sm 2 6 chz < 0.041 , (1) 

where in the normal hierarchical scheme, mi << to 2 << 
7713, we have Am 2 ol = Am^, Am 2 tm = Am^ and so i — 

#12, Oatm — $23 7 "chz = $13- 

In order to explain the above neutrino mass differences 
without introducing extremely small and unnatural neu- 
trino Dirac Yukawa couplings, the general practice is to 
add superheavy right handed Majorana neutrino fields, 
Ni, and rely on the seesaw mechanism for generating sub- 
eV light neutrinos: 

C = C SM + C Y (N, L, 0) + M%NiNj/2 , (2) 

where, in the minimal setup, Cy(N,L,c/)) — Y^LiHNj, 
with Li — {vjLj-jL) and H is the only Higgs doublet 
that couples to neutrinos. Then, the light neutrino mass 
matrix is given by the see-saw mechanism as: 

m„ = -m D M N l m T D , (3) 

with mrj = vY v and v =< H >. 

In this work we wish to propose a 2HDM, which we will 
name the "2HDM for the 3rd generation" (3g2HDM). 
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The model relates the neutrino sector to the up-quark 
sector, and by giving the third generation right-handed 
neutrino similar status as the top quark it provides a nat- 
ural framework for accommodating the observed pattern 
of neutrino masses and mixing angles. In addition, our 
3g2HDM assumes a specific structure in the Yukawa sec- 
tor which should be viewed as an effective low-energy 
parametrization of the underlying short distance the- 
ory. This concrete phenomenological framework for the 
Yukawa sector allows us to successfully address both the 
heaviness of the top-quark and the apparent hierarchical 
structure in the neutrino sector. 

In particular, we extend the idea of the so called 
"2HDM for the top-quark" (t2HDM) j to the leptonic 
sector. In the t2HDM one assumes that </>2 [the sec- 
ond Higgs doublet with a much larger vacuum expec- 
tation value (VEV)] couples only to the top-quark, while 
the other Higgs doublet <f>i (with a much smaller VEV) 
couples to all the other fermions. The large mass hi- 
erarchy between the top quark and the other quarks is 
then viewed as a consequence of V2/V1 = tan/? >> O(l), 
which, therefore, becomes the "working assumption" of 
the t2HDM. 

The t2HDM has several other notable features, such as 
enhanced H + bc and H°cc couplings, new flavor changing 
H°tc and H°tu interactions and new CP phases. Thus, 
the t2HDM can influence CP- violation in B physics 
flavor changing Z-decays and the production of a 
Higgs in association with a c-jet or a b-jet in hadron 
colliders 0. 

Following the idea of the t2HDM, we propose the fol- 
lowing Yukawa interactions in the 3g2HDM: 



C Y = - Y e L L ^B-Y d QL<i>id R 

- Y'LLfaN -Y»L L 4> 2 N + h.c. 



(4) 



where TV are right-handed Majorana neutrinos, Q and L 
are the usual quark and lepton doublets and 



yU 










c u < L 
,0 c u ' L 



,(5) 



such that, in both the quark and leptonic sectors, </>2 
couples only to the third generation right-handed up- 
fcrmions. The Yukawa texture in (JSJ yields the following 
Dirac neutrino mass matrix (dropping the superscript u): 



niD 




a b ~ 
a b ctp 
5b ctp , 



(6) 



where tp = tan/3 = v^jvx. 

Note that we have assumed the texture zeros 
(Y™' v )i3 = and (Y£' v )ai = 0, which gives a Dirac neu- 
trino mass matrix in which the upper right and lower 
left entries vanish. This is required in order to simplify 



our discussion below and to avoid fine tuning in the neu- 
trino Yukawa sector when confronted with the experi- 
mental limits on neutrino masses and mixings (in fact, 
it is sufficient to set only (Y 2 l ')i3 — > in order to avoid 
fine tuning). The Dirac mass texture in © is similar to 
the one which arises in the minimal approach suggested 
m 0- j that is, from the requirement of the most eco- 
nomic addition to the SM that can accommodate both 
neutrino oscillation data and the observed baryon asym- 
metry in the universe. Clearly, for such texture zeros 
to be natural, they have to be protected by underlying 
symmetries of the short distance theory. For example, 
SO(10) GUT theories dictate a symmetry between the 
up-quark and the Dirac neutrino mass matrices which, 
along with the empirically known properties of the up- 
quark mass matrix, yields the texture zeros in ©, see 
e.g., 0. In fact, as was shown in |10|. there always exists 
a scalar sector such that any texture zero in the fermion 
mass matrices can be enforced by means of an abelian 
symmetry once the scalar sector is enlarged to the extent 
required. Another alternative and attractive framework 
for the zero entries in Y\ and Y2 was suggested in 8] : such 
texture zeros may easily result from a vanishingly small 
overlap between the wave functions of N and L with the 
corresponding indexes in extra-dimensional models. The 
above several possible ways to accommodate the desired 
texture zeros [in 0] do not give us any further insight 
regarding the underlying reason for the structure of the 
Yukawa couplings in our model. Nonetheless, they do 
motivate us to study the implications and phenomenol- 
ogy of the 3g2HDM with the texture zeros assumed in 

We note also that, without loss of generality, further 
simplifications were taken in that allow us to give a 
more transparent analytical derivation of the neutrino os- 
cillations pattern in our 3g2HDM. In particular, we have 
chosen similar Yukawa couplings within each generation: 
(Y?' v J a (i = 1,2) ~ a«>", (Y?>") i2 (i = 1,2,3) ~ b u > v and 
{Y 2 w ' v )a{i — 2,3) ~ c u ' v . In addition, a real parameter 
of O(l) (the parameter 5) was added to the (Yi)23 entry 
in JSJ . We emphasize that these simplifying assumptions 
have no effect on our main results below, aside from sim- 
plifying the analytical formulae given below. 

In the basis where Mm is diagonal (one can always 
choose a basis for the fields N such that Mm is diagonal) , 
Mm — M ■ diag(eMii e M2, (M3), we then obtain from the 
seesaw mechanism in Q: 



(7) 




where 



and 



a 



b 2 

ZM2 



(-Ml 



C 2 t 2 
CM3 



(8) 



(9) 
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In the limit uj >> e, e, 8, the physical neutrino masses 
and mixing angles arc then given by 1 1] : 



tan 2#23 

771l 



2ruj 



-, tan 20i2 



2- 

/' 



e(6 2 

eml { 1 - g sin 20 12 + / sin 2 12 } 
to 2 — em° { 1 + g sin 20 12 + f cos 2 6 12 } , 
777,3 ~ 2ujni Q u , 
where we have defined the parameters: 

S\ . S 2 - 26 - 
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e(S- 



2 3 / 2 



V2r 



f 



2r 



(10) 



(11) 



In the following, we will assume normal hierarchy for the 
light neutrinos: 7773 >> m 2 >> mi, therefore leading to 

. Thus, setting for 



m 3 



V Am lt m and m 2 ~ 



example 2mP v ~ ^/ Am 2 o; , it follows from the expression 



for 7773 in HI L)| l that uj 



V Am2 at m /Amj ol . 



Note that, using the definitions for w in JS}, for raJJ in 
(O and the fact that 7774 ~ ctp x (vi/^/2), our model 
yields the following seesaw triple-relation between the 
heaviest light neutrino, the heaviest right-handed Ma- 
jorana neutrino and the top quark: 



7773 



(12) 



Moreover, from © with 2m" ~ \J Am 2 o; we obtain the 
typical mass scale (M) of the heavy right handed neutri- 
nos: 



M 



vol 



10" GeV , (13) 



'sol 



where, in the last equality, we have taken mt ~ v 2 j\[2 
and tp ~ O(10), which are the "working assumption" 
values within the 3g2HDM. 

In the following numerical analysis, we will set uj 



V Am i tm /Ami ol = 5.18 and 



= 4.53 • 



v wviivi sui ** v but, 1 

10~ 3 eV, which correspond to the best fitted values [2|: 
Am 2 tm = 2.2 • 10~ 3 [eV] 2 and Am 2 sol = 8.2 • lO" 5 [eV] 2 . 
We note, though, that the best fitted set of input param- 
eters (obtained by performing a minimum \ 2 fit of our 
model to the experimentally measured values of atmo- 
spheric and solar neutrino masses and mixing angles Q) 

is: pPf u ~ 5.34, e - 0.57, S 1.28 and r ~ 1. 

In Fig. we give a scatter plot of the allowed ranges 
in the #13 — #23 and #13 — 0\ 2 planes, i.e., subject to 
the 3cr limits listed in Q). This is done by randomly 
varying the three input parameters e, r and 5 (requiring 
uj >> e, e, S) with a sample of 3 • 10 6 points. We see 
that the 3g2HDM predicts #13 to lie within (in radians) 



with the best fitted value at #13 ~ —0.011 and #23 ~ 
43.9°. Notice that this is a rather restrictive prediction 
for the atmospheric mixing angle since the experimen- 
tally allowed 3cr range for #23 spans over about 20°, i.e., 



35° < e 2 x / < 55° [see ©] 



2 
1 

M 

5 -1 

rH , 




• * 1 . » . 



41 42 43 44 45 46 

e 23 [deg] 




27 29 31 33 35 37 39 



e i2 [deg] 



FIG. 1: The allowed ranges in the #13 —823 (left) and 613 — 6*12 
(right) planes, for to ~ 5.18 and ml ~ 4.53 TO" 3 eV. 5ee a/so 
text 

To better understand what is enforcing maximal atmo- 
spheric mixing and a very small #13 (as is evident from 
Fig- EI we note that from 1)10(1 we have: 



, e{5 + r) „_ x 

?13 ^ X UJ 



V2 



723 



1 

— tan 
2 



e(<5 2 - r) 



where 



2ru> ~ 2r 4 



Am 



» 1 



,(15) 



(16) 



since by definition r > 1. Therefore, for values of e, (5 
and r of 0(1), we have: 



(17) 



as is typical to grand unified (or quark-lepton unified) 
theories [l-l . Furthermore, from (|15|) we have (taking 
<5 < — 1, see discussion below and Fig. EJ) : 




723 




'13 



(18) 



as can be seen in Fig. ^ 

In Fig. |21 we give a scatter plot of the allowed regions 
in the S ~ e plane, in two limiting cases |l2j| : 



-0.05 £ 9 13 £ 0.035 and the atmospheric mixing angle case I: r ~ 1 (e ~ e), corresponding to -f^ « 
#23 to be at most a few degrees away from maximal. In 
fact, a minimum x 2 analysis with respect to O13 and 



'23 



yields: 



0.017 < 0i3 
42.9° < 23 



< 0.021 

< 45.2° 



99% CL , 
99% CL , 



(14) 



case II: r ~ 2 (e — > 2e), corresponding to -2— ~ . 

V " ^ to £M1 EA72 

We see that case I (corresponding to the best fitted value 
for r 12}) is compatible with neutrino oscillation data 

1. while the allowed 
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case I: r=l 



case II: r=2 



OJ 0.6 




Furthermore, we find that our Ansatz (|19fl yields the fol- 
lowing masses for N 2 and Ni [for r ~ 1 and for e ~ 0(1) 
in accord with the constrains coming from neutrino os- 
cillations, see Fig. |2]: 



M N2 ~ l(T 2 Af , M Nl » 10~ 6 M 



(21) 
10 11 



FIG. 2: The allowed ranges in the e — 5 plane, for cases I 
(left) and II (right). See also caption to Fig. 



TABLE I: Mixing angles and neutrino masses in cases I and 
II with e = 0.5 and e = 0.75, respectively. For both cases 
5.18 and m° ~ 4.53 ■ 10~ 3 eV. 



S = —1.5, oj 



which, for M ~ 10 13 GeV [see JT3J], gives M N . 
GeV and M Nl » 10 7 GeV. 

Let us now examine how the 3g2HDM with our quark- 
-2.5 -2.3 -2.1 -1.9 -1.7 -1.5 -1.3 lepton similarity Ansatz in l|19l) and with r ~ 1 and 

e ~ 0(1) fits into the mechanism of leptogenesis [l6| . 
Leptogenesis is an additional attractive feature, associ- 
ated with the seesaw mechanism, since it allows the pos- 
sibility of generating the observed baryon excess in the 
universe from a lepton asymmetry, driven by decays of 
the heavy Majorana neutrinos [l4. ITEL H^|. In particu- 
lar, a CP-asymmetry, ejf f , in the decay Ni — ► £<j>j can 
generate the lepton asymmetry (see e.g., fl4ll5| ): 





S = -1.5, uj ~ 5.18, m°„ ~ 4.53 ■ 10" 3 eV 


case (e) 


#23 


012 


013 


mi [eV] 


m2 [eV] 


m 3 [eV] 


I (0.5) 
II (0.75) 


43 u 
44.7° 


29.5 U 
35.9° 


-1.14 u 
0.69° 




3.5 ■ 10~ 4 


0.0093 
0.0092 


0.047 
0.048 



range for case II is 0.5 £ e £ 0.8 with -2.4 < <5 < - 1.4. 
As an example, in table [I] we give the neutrino mixing 
angles and masses in cases I and II, for some specific 
values of the allowed parameter space. 

Let us now examine the consequences of our model 
on the mass spectrum of the heavy Majorana neutrinos, 
taking into account the above constraints coming from 
neutrino oscillation. For dcfinitcncss, we will consider 
case I (i.e., r ~ 1, which we find as the best fitted value 
with respect to oscillation data), and based on quark- 
lcpton similarity (perhaps motivated by GUT scenarios, 
see e.g., [t|), we will take the following Ansatz for the 
neutrino Dirac Yukawa couplings: 



O(10" 3 ), 6~ 0(10^), c~ 0(1) 



(19) 



Indeed, this choice follows from the assumption that the 
Dirac mass matrices of the neutrinos and up-quarks ex- 
hibit a similar hierarchical structure. This Ansatz is, 
therefore, natural in the 3g2HDM under consideration, 
since this model is constructed in order to simultaneously 
explain the large top mass and the hierarchical struc- 
ture of the neutrino masses, by relating the quark sector 
to the leptonic sector. That is, taking V\ ~ O(10) and 
tp ~ O(10), this Ansatz follows from the up-quark sector; 
a v v% ~ 0(m u ), b u vi ~ 0(m c ) and m t ~ O(c u vita,n0), 
if one assumes the quark-lepton similarity: a" ~ a u ~ a, 
b u ~b u - 6 and c u ~ c v - c. 

It is then easy to extract the mass spectrum of the 
heavy Majorana neutrinos in our Ansatz, subject to the 
constraints given in Q). In particular, from l|12|) we have: 



M 



N 3 



10 15 GeV 



(20) 



n L /s = e Ni Y Ni (T » M Ni )rj 



(22) 



where = n^Js, njsf i being the number density 

of Nj and s the entropy, with Yat 4 (T >> MjvJ = 
135C(3)/(4-7r 4 g*) and g* being the effective number of 
spin-degrees of freedom in thermal equilibrium. Also, rj 
is the "washout" parameter (efficiency factor) that mea- 
sures the amount of deviation from the out-of-equilibrium 
condition at the time of the N decay. 

The lepton asymmetry n^/s can then be converted into 
a baryon asymmetry through nonperturbative spheleron 
processes. The conversion factor is n B/s = 

-{8N G + AN H )/(22N G + 13 N H ) x n L /s] where N G and 
Nh are the number of generations and Higgs doublets, 
respectively. Thus, in our case, i.e., N G — 3, Nh = 2 
and ^ 100, we obtain: 



n B /s~ -1.4x10 3 e Ni r) 



(23) 



As seen from (|2UJ) and l(2T]l. our 3g2HDM can lead to a 
hierarchical mass spectrum for the heavy Majorana neu- 
trinos, MjVj << Mpf 2 << Mjsr 3 , within a large portion 
of the allowed parameter space. In particular, imposing 
M Nl « M N2 , from and with M - 10 13 GeV, we 
get: M Nl - 10 8 - 10 10 GeV. With such an hierarchical 
mass spectrum, only the CP-asymmetry produced by the 
decay of N± survives. 

In our model N\ can only decay to the "light" doublet 
4>i, as can be seen from (0J and (J5J. Therefore, the CP- 
asymmetry £]\[ 1 is similar to the one obtained in single- 
Higgs models in the limit Mjv x << Mn 2 , Mry 3 (see e.g., 

my- 



ejVi — 



16?r ^ 

i=2.3 



Im 




M Nl 


( 




M N . 



(24) 



where Y" is the Yukawa coupling of N\ to cf>i given in 
(JSJ). Note that, since N 3 does not couple to <f>x in our 



5 



model, it does not contribute to the sum in l|24[l above. 
Therefore, is generated only from the i — 2 term in 
J2IJ, giving: 



36 2 M Nl 
"8n M 



N 2 



sin2(6» & 



(25) 



where the CP-phases arise from the possible complex en- 
tries in Yj": a = \a\e w " and b = \b\e ie » . Thus, using O 
and e ~ b 2 /e M2 (as implied by a fit of the 3g2HDM to 
neutrino oscillation data, i.e., r ~ 1), we obtain: 



ti 



16tt 



eMjv lS in2(0 6 -0 a ) . (26) 



It is interesting to note that the above CP-asymmetry 
is a factor of ~ £| larger than the CP-asymmetry ob- 
tained in models with one Higgs doublet or for that mat- 
ter in SUSY-like models where the condensate of a single 
Higgs is responsible for generating the Dirac neutrino 
mass term. Thus, in our 3g2HDM the enhancement fac- 
tor is of 0(100) since tp ~ O(10). To see that we can 
rewrite cn x in Ij25(l as follows: 



sin2((9 fc -6 a 



2e , 
x — ti 



(27) 



where in our model e ~ 0.5 and 



5 are fixed by 



oscillation data and £ 7 ^ x is the maximum of the CP- 
asymmetry in models with one Higgs doublet (see e.g., 



asyi 

m 



)■■ 



3 M Nl m atm 

167T V 2 



(28) 



with the atmospheric neutrino mass m a 



in (IKjp and v = \]~v\ 



= m% given 
246 GeV. As mentioned 
above, this enhancement is possible due to the presence 
of the second Higgs doublet, since in this case the light 
neutrino mass spectrum is determined by v\ (the VEV 
of the lighter Higgs field) and not by the EW-scale VEV 
v (as in the usual scenario with one Higgs doublet), and 
the hierarchy m a tm >> m-soi is dictated by the large 
tan/3. In this way, the CP-asymmetry no longer suffers 
from a direct relation to the Dirac neutrino mass term 
and is, therefore, not directly constraint by oscillation 
data. This observation was also made by Fukuyama and 
Okada in 0] wn0 have investigated leptogenesis within 
a generic two Higgs doublet model. As will be shown 
below, this enhancement of the CP-asymmetry allows us 
to reproduce the observed baryon asymmetry in the uni- 
verse with Mjvj as low as 10 9 GeV, which in turn re- 
laxes the lower bound on the reheating temperature in 
the early universe to T RH > 2 ■ 10 8 GeV, [3. 

The efficiency factor, r], which measures the amount 
of lepton asymmetry left from the CP-asymmetry , 
is calculated by solving the appropriate Boltzmann 
equations. The key parameter entering the Boltz- 
mann equations is the "decay parameter" , defined 



as [13: K = T Nl /H(T ~ M Nl ), where T Nl = 
(Y{ t Yf) 11 M Nl /8ir is the total decay width of Ni and 
H(T) = y/i-n-Zg* /45T 2 /M P i anck is the Hubble expansion 
rate. In particular, using (JSJ, the decay parameter K is 
given by (for ~ 100): 



K 



x 10"V 



Mpi 



anck 



M 



(29) 



which, for a — 10~ 3 [i.e., Ansatz 0], gives 6 £ K £ 60 
if e.g., M Nl - 10 9 - 10 10 GeV. In this range of values for 
the decay parameter (corresponding to the "mildly strong 
wash out" regime), the relation rj ~ 0.5/ K 1 - 2 constitutes 
a good fit to the numerical solution of the Boltzmann 
equations, see e.g., P. Di Bari in |l5j . Thus, using the 
above fit for rj and the CP-asymmetry in l|26() , the baryon 
asymmetry in (|23|) becomes: 



n B 



10" 



17 "P 



tU Ami 



2m 2 



sm2(0 b -9 a ){3O) 



This has to be compared with the observed baryon to 
photon number ratio ns/n-y ~ 6 x 10~ 10 19], implying 
Ub/s ~ 8.5 x 10~ n . For example, taking e ~ 0.5 (see 
Fig.0 and Am 2 sol - 8.2 • 10~ 5 eV 2 , along with tp - 10 
and mt ~ 170 GeV, (|3H|l reproduces the observed baryon 
asymmetry for e.g., M^ x ~ 10 10 GeV and sin 2(9b — a ) ~ 
0.1, or for MjVj ~ 10 9 if CP is maximally violated in the 
sense that sin 2(06 — 6 a ) ~ 1. 

Summarizing, we have presented here a simple exten- 
sion to the Standard Model, in a grounds-up approach, 
which leads us to a framework for neutrino masses and 
mixings with considerable predictive power. This can 
clearly have significant impact on on-going as well as 
planned experiments. To briefly recapitulate, we have 
constructed a 2HDM in which the second doublet, with 
a much larger VEV, couples only to the third genera- 
tion right-handed up-fermions (the top-quark and the 3rd 
generation right-handed neutrino) and the other doublet 
couples to all other fermions. Thus, this model is a possi- 
ble effective low energy parametrization of an underlying 
short distance theory which envisions a close relation be- 
tween quark dynamics and neutrino physics. The key 
parameter of this 2HDM is tan /? which is assumed to be 
of O(10) in order to naturally accommodate a large mass 
for the top quark. We have shown that the large value 
of tan P in this model is directly responsible for success- 
fully reproducing the observed neutrino oscillation data, 
predicting a very small 9 13 : -0.017 < 9 13 < 0.021 at 99% 
CL, a very restricted allowed range for the atmospheric 
mixing angle: 43° Z 9 23 £ 45° at 99% CL, as well as suc- 
cessfully reproducing the observed baryon asymmetry of 
the universe through leptogenesis. In particular, in our 
model, the CP-asymmetry in the N\ decays that drives 
leptogenesis is larger by a factor of 0(tan 2 (3) - O(100) 
than the CP-asymmetry obtained in models for lepto- 
genesis with one Higgs doublet or in the MSSM. This 
enhancement allows to relax the lower bound on 
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and accordingly also the lower bound on the reheating 
temperature of the early universe. 
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